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New results concerning the production of neutral strange particles, K^ and A 
in Ni + Ni collisions at 1.93^1 GeV, measured with the FOPI detector at GSI- 
Darmstadt, are presented. Rapidity density distributions and Boltzmann slope pa- 
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rameter distributions are measured in nearly the full phase space of the reaction. The 
observables are compared to existing and proton data. While the data agree 
with previously reported measurements, the A distributions show a different be- 
havior relative to that of protons. The strangeness balance and the production yield 
per participating nucleon as a function of the centrality of the reaction are discussed, 
for the first time at SIS energies. 

PACS numbers: 25.75.-q, 25.75.Dw 



I. INTRODUCTION 

Strangeness production and propagation in relativistic heavy-ion collisions are active 
research topics for experimental and theoretical nuclear physics since they are expected to 



ows to 



provide interesting opportunities for studying hot and dense nuclear matter. It al 
address fundamental aspects of nuclear physics such as the nuclear Equation of State 
and the question whether hadronic properties undergo modifications in such an environment 
0, y, ^. Moreover, it is also of crucial interest to improve our understanding of the reaction 
mechanisms governing these collisions f^]. The field of heavy-ion physics is also of great 
importance for astrophysics, in particular to investigate the characteristics of the core of 
neutron stars Q]. 

The SIS energy range GeV is best suited to study the in-medium properties of 

strange particles since they are produced below or close to threshold. The density of the 
nuclear system created during the reaction is expected to reach up to three times normal 
nuclear matter density and its temperature is of about 90 MeV Q . Theoretical works 
predict that chiral symmetry can be partially restored under those conditions, leading to 
changes of hadron properties [ll] affecting both, production and propagation. Indications 



for in-medium modifications of charged 



^aons have been already observed experimentally 



with data from FOPI and KaoS 3. 

Neutral strange particles such as and A are of great interest because of their associated 
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production and the fact that they are not influenced by the Coulomb interaction. The study 
of their properties may probe the in-medium potential, which is expected to be weakly 
repulsive for (K^) uM attractive for A [16]. It is worth to point out that such 
opportunity cannot be offered by studies of hypernuclei [17] which allow to investigate the 
A-nucleon potential at normal nuclear matter density. 



Detailed studies of the properties of these neutral strange particles have been per 



at t 
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he LBL BEVALAC and the BNL Alternating Gradient Synchrotron (AGS) (if 



brmed 



22| . In this paper new results on the production of and A in Ni + Ni collisions 



at 1.93^4 GeV, near threshold (which is 1.58 GeV for free A^A^ collisions), are presented. 
Differential production yields and inverse slope parameters from close to the full phase 
space are discussed. In particular the high statistics collected during the experiment allows 
for the first time to investigate the centrality dependence of the and A production yields 
at SIS energies. 

The paper is structured as follows. Section [TTl consists out of a short description of 
the apparatus. Section [TTIl is devoted to the event characterization, namely centrality and 
azimuthal orientation of the reaction plane. The reconstruction method is detailed in section 
HVl Section [Vl addresses the experimental results, yields and effective (Boltzmann) slope 
parameters of and A, which are also compared to existing data of and proton. 
Subsequently the strangeness balance and the centrality dependence of and A production 
are discussed. Finally, a summary and an outlook are given in section IVll 



II. EXPERIMENT 



The experiment was performed at the GSI Schwerionen Synchrotron (SIS) in Darmstadt 
(Germany) by bombarding a 28^1 beam of 1.93^1 GeV on an enriched target (>95%) of 28^1. 
The target thickness was about 360 mg/cm^ corresponding to an interaction probability of 
1.5% and the average beam intensity was 4 - 5x10^ ions/s. 

The FOPI setup is an azimuthally symmetric apparatus made of several sub-detectors 
which provide charge and mass determination over nearly the full An solid angle. The 
central part (23° < 9iab < 113°) is placed in a super-conducting solenoid and comprises a 
drift chamber (CDC) surrounded by a barrel of plastic scintillators for polar angles between 
32° and 72°. The mass of particles measured in the CDC is determined using magnetic 
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rigidity and energy loss. The forward part is composed of a wall of plastic scintillators 
(1° < 9iab < 25°) and another drift chamber (Helitron) mounted inside the super-conducting 
solenoid. The plastic wall is divided in two parts: the inner plastic wall and the outer plastic 
wall (PLAWA) which cover polar angles 1° < 9iab < 6.6° and 6° < 9iab < 25°, respectively. 
The plastic wall provides charge identification of the reaction products, combining time 
of fiight and specific energy loss informations. For the present analysis, the PLAWA and 
the CDC were used. More details on the configuration and performances of the diflFerent 
components of the FOPI apparatus can be found in 

During the experiment, more than 110x10^ events were recorded under "central" and 
"medium-central" triggers, corresponding to different thresholds on the multiplicity of 
charged fragments measured in the PLAWA. Most of the events (100x10^) were registered 
under the "central" trigger which accepts events covering up to about 20% of the total 
reaction cross section. 

In previous experiments Q, Q, Q, 



l24l | the target was located at a diflFerent position 
+40 cm in beam direction which led to a better acceptance around target rapidity. With 
our present target setup, the whole backward hemisphere is covered. 

III. EVENT CHARACTERIZATION 
A. Centrality Selection 

The events are sorted according to their degree of centrality by imposing cut conditions 
on the total charged particle multiplicity measured in both, PLAWA and CDC (i.e. MUL). 

Four event classes are used for the investigation of the centrality dependence of the 
and A production yields (Tab. H labeled "MULl" to "MUL4"). The characteristics of the 
central event class used to study the and A diflFerential distributions are provided in the 
last row of Tab. [J (labeled "Central"). The corresponding ranges (from (7geo,min to (7geo,max) 
and the portions of the geometrical cross section of the centrality classes (Aageo) are listed. 
The mean geometrical impact parameter (fe^eo) and its r.m.s. deviations are determined 



assuming a sharp-cut-off approximation [25|]. We calculate the average number of nucleons 



in the fireball (Apart) using the geometrical model described in [26]. 
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^ CjCO^TTLOLX 




(bgeo) 




Class 


(mb) 


(mb) 


(mb) 


(fm) 




MULl 





98 


98 


1.7 ± 0.4 


94 ± 6 


MUL2 


98 


389 


291 


3.6 ± 0.4 


63 ± 6 


MUL3 


398 


679 


290 


4.7 ± 0.3 


43 ± 5 


MUL4 


679 


1005 


326 


5.5 ± 0.3 


32 ± 5 


Central 





350 


350 


3.1 ± 0.4 


71 ± 6 



TABLE I: Cross section, mean geometrical impact parameter, and number of participants for 
Ni + Ni collisions at 1.93A GeV for five different centrality selections (see text for details). 

B. Reaction Plane 



When applying the event-mixing technique for the combinatorial background determi- 
nation one needs to rotate all events into the reaction plane (see section IV). This plane 
is estimated on an event-by-event basis according to the standard transverse momentum 
procedure devised in [27] which allows to construct the event-plane vector 

Q ^UJ^Uiy . (1) 

The sum runs over all charged particles measured in the PLAWA and the CDC. The unit 
vector Ujy is parallel to the particle transverse momentum (i.e. Uj^ = (cos^^^^, sin Lpiy)^ where 
(fiy is the particle azimuth), and uUiy is a weight to improve the resolution, depending on the 
center-of-mass (cm.) rapidity yc.m. = Viab — 0.5 • yp (the subscript p refers to the projectile): 
uUiy = —1 for ycm. < —0.2, cj^ = +1 for ^c.m. > +0.2 and uji^ = otherwise. The azimuth 
of Q is an estimate of the azimuth of the reaction plane. Typically, the resolution on the 



reaction plane determination estimated as in 



271] from randomly chosen sub-events, is about 



30° for central events {(Jgeo ^ 350 mb, bg^o < 4.0 fm). 



IV. PARTICLE RECONSTRUCTION 



Long-lived neutral strange particles like the or the A are identified via their charged 
decay particles in the CDC. and are produced in equal amounts due to the nature of 
the weak interaction (neglecting CP violation). Because of its lifetime and decay channels 
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the cannot be measured with the FOPI apparatus, whereas the can be reconstructed 
from its decay into tt" and tt^ (branching ratio = 68.6%). Similarly, the reconstruction of 
the A is done from the final state containing tt" and p (branching ratio = 63.9%). A possible 
contribution to the observed A comes from the TP which decays exclusively into a A and a 
photon. Since the photon cannot be measured with the FOPI detector, A and TP cannot be 
disentangled and are referred to as A in the following. 

The long lifetimes of the weakly decaying [cr = 2.68 cm) and the A (cr = 7.89 cm) 
causes a sizable fraction of these particles decaying somewhat away from the primary vertex. 
The precision of the track reconstruction in the CDC is sufficient to resolve these secondary 
vertices. Therefore the signal-to-background ratio is improved significantly by rejecting 
particles coming from the primary vertex. 

A. Reconstruction Method 

The search for and A starts with the selection of possible final state particles. Three 
cuts for each track in the CDC are applied. The first one rejects particles with transverse 
momenta lower than 80 — 100 MeV/c which are spiraling in the CDC. The particle species is 
selected using the correlation of the particle mass, which is determined from the measured 
specific energy loss, and the total momentum. The distance of closest approach between 
track and primary event vertex in the transverse plane is the third quantity. It is used to 
enhance the fraction of particles coming from secondary vertices in the data sample. 

For a jet-type drift chamber like the CDC it is advantageous to pursue a reconstruction 
strategy which first focuses on the transverse plane (x — y) in order to find intersection points 
of possible decay particle tracks as illustrated in the upper plot of Fig. [H Once such a point 
is found, further cut quantities are determined. The transverse distance to the primary 
vertex rt and the closest approach of the back-extrapolated mother particle fiight path do 
are used to suppress the background. Then the quality of the signal can be improved by 
taking into account the longitudinal distance Az of the decay particle tracks as shown in 
the lower plot of Fig. [H Afterwards, the invariant mass of the particle pair is calculated 
from its four-momenta at the intersection point. 
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B. Invariant Mass Spectra and Phase-Space Distributions 

The invariant mass spectra of 7r~7r+ and 7T~p pairs for a centrahty selection of the most 
central 350 mb of the geometrical cross section are shown in Fig. [2l A total of about 
47x10^ events were analyzed. The combinatorial background was modeled with the event- 



mixing method (28[]. For that purpose the two decay particle candidates are taken from 
two different events with the same centrality. In addition, the two events are projected 
into the same reaction plane. Afterwards, a normalization is applied to the mixed-event 
background in order to fit the combinatorics in the invariant mass range 0.56 — 0.80 GeV/c^ 
and 1.15 — 1.25 GeV/c^ for and A, respectively. The shape of the resulting mixed-event 
background fits nicely to that of the combinatorial background and is indicated by the shaded 
area in Fig. [2l After background subtraction, the remaining peaks in the mass spectra are 
fitted with Gaussians. Within an interval of ±2a around the center of the peak, about 33000 
and 61000 A candidates are found for the applied set of cuts in the analysis. The signal- 
to-background ratios are of 0.9 (K^) and 2.1 (A). The (Gaussian) widths of 15.8 MeV/c^ 
(K^) and 4.1 MeV/c^ (A) are resulting from the detector resolution. 

The background-subtracted phase-space distributions of those particle candidates are 
presented in Fig. [3] where the transverse momentum pt is plotted versus the center-of-mass 
rapidity yc.m.- The solid lines denote the limits of the CDC's polar angular acceptance 
of 23° and 113° in the forward and backward directions, respectively. Both distributions 
illustrate the nearly complete coverage of the backward hemisphere. By using the symmetry 
of the colliding system, most of the phase space becomes accessible. In addition, the 
distribution shows a fraction of candidates even coming from outside the CDC acceptance. 
This is possible because both decay particles can still reach the CDC acceptance due to the 
high excess energy in the decay process. 



C. Acceptance and Efficiency Correction 

In order to extract global quantities like the particle yields from the data, the losses 
due to acceptance and efficiency have to be corrected. For a complex apparatus like the 
FOPI detector these corrections have to be determined by an extensive GEANT simulation 
modeling all detector components. The particle distributions of heavy-ion collisions are 



generated by the IQMD model [9|, l29|]. The neutral strange particles (X^, A) were added 
separately using a momentum distribution given by the Siemens-Rasmussen formula jsol 
which describes a system expanding with a temperature T and a radial expansion velocity 
f3. The choice of the parameters (T = 90 MeV^ (3 = 0.3) is suggested by previous results 
obtained for tt", protons and deuterons |lol]. and A particles were then embedded into the 
IQMD output with one and one A per event. A full simulation of the detector including 
resolutions in energy deposition and spatial position, front-end electronics processing and 
hit tracking, is performed and the resulting output is subject to the same reconstruction 
procedure as the experimental data. Spectra of all relevant cut quantities were compared 



between simulated and experimental data 
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321]. Since no significant diflFerences were 



found, one can conclude that the detector is properly described by the simulation. 

Finally, the reconstruction eflSciency is determined by computing the ratio of the number 
of reconstructed particles in the simulation to the number of particles embedded into the 
input events. This is done in bins of rapidity ^c.m. and transverse mass — mo, where 

= ^/'^l + Pt mo is the rest mass of the particle. Thus, the correction method is 
independent of the choice of the source temperature for the embedded strange particles. 

The reconstruction efficiencies for (open circles) and A (full squares) are plotted in 
Fig. m The numbers are not corrected for the branching ratio of the respective decay chan- 
nels. Due to the fact that a minimum transverse momentum of 90-100 MeV/c is requested 
for the decay particles, the efficiency is low for small values of rrit — mo (< 100 MeV/c^). At 
higher transverse masses the efficiencies are approximately independent of i/cm. and rrit — mo 
with maximum efficiencies of about 7% for the and 5% for the A, respectively. In order 
to apply a smooth efficiency correction to the data, which is less sensitive to ffuctuations in 
the high-mt range, the following function was fitted to the reconstruction efficiencies: 

f = a ' tanh(6 • {rrit ^o)) + c • {rrit '^o)^ , (2) 

where a, b and c are free fit parameters and d has a value of 0.45 and 0.2 for and A, 
respectively. This function was then scaled by the branching ratio of the decay channel and 
applied to the measured transverse mass spectra. 
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V. RESULTS 

In this section, the transverse mass spectra of and A after background subtraction 
and efficiency correction, are fitted with a Boltzmann-type distribution in order to extract 
the rapidity dependencies of the Boltzmann temperature, the rapidity density distributions 
and ultimately, the effective temperature and the total production yield in central (350 mb) 
Ni + Ni reactions at 1.93A GeV. The distributions obtained for and A are compared to 
existing data on and protons, respectively. The large statistics of reconstructed and 



A allows for the first time to discuss the total production yield as a 



"unction of the centrality 



of the reaction, as it has been done at AGS energies [l8|, ll9|, l20l. l22[|. 

In the following, the results include only statistical errors. Note that the differential yield 
distribution of the was scaled by a factor two in order to account for the unmeasured 
contribution from the and is therefore referred to as K^. 



A. Transverse Mass Spectra 

Figure [5] shows the transverse mass spectra rrit — mo of and A particles, after back- 
ground subtraction and correction for efficiency, branching ratios and the unmeasured 
contribution, for several rapidity bins indicated in the figure. The rapidity range 
— 1.0 < yc.m. < 0.0 is covered, hence the full phase space is accessible by exploiting 
the symmetry of the Ni + Ni system. In order to compile all spectra in one plot they are 
multiplied by 10^ starting from the lowermost spectrum (n = 0) to the uppermost spectrum 
(n = 9). 

If particles are emitted from a thermal Boltzmann-like source, their behavior is then 
described by the following function: 

i, ,, =A.e.p -^'^l-'^'^ (3) 

mf d{mt - mojdyc.m. Tb 

within a narrow window of rapidity dyc.m.- Both, the inverse slope parameter [Tb] and the 
integration constant {A) depend on the rapidity. This allows to extract the rapidity distri- 
butions of the inverse slope parameter Tb and the rapidity density distributions dN/ dyc.m. 
by integrating the fitted function from — mo = to oc. In the logarithmic representation, 
our measured spectra exhibit a linear decrease with increasing rrit — mo and are described 
reasonably well by Eq. [3l 
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B. Inverse Slope Parameters 

The inverse slope parameter Tb yields information on the particle temperature at freeze- 
out for a thermalized system. The Boltzmann temperature Tb has a simple dependence on 
the rapidity for an isotropically emitting source: 

Tb= If y (4) 

where Tgjj = TB{yc.m. = 0) is the effective temperature, i.e. the slope parameter at mid- 
rapidity. 

The rapidity dependencies of Tb are depicted in Fig.[6]for various particle species. The full 
circles and squares indicate measured data points for and A while open symbols denote 
points reflected with respect to mid-rapidity. The star symbols correspond to previously 



measured distributions of from [2J] and to protons from this experiment. For the 
an effective temperature of 114 ± 1 MeV is extracted from the flt using Eq. [4] which 
agrees very well with the results from [421]. The Tb distribution of is found to be 



slightly higher than the one measured for K^^ within the overlap region. Note that 
were measured directly, i.e. via time-of-flight and momentum, and that in previous FOPI 
experiments the target was in its nominal position which explains the different rapidity range 
for K^. Comparing the combined KaoS and FOPI results for the to the suggests 
that they have very similar freeze-out conditions. 

The lower part of the Fig. [6] displays the corresponding results for the A with an the 
effective temperature of 119 ± 1 MeV. The A have a somewhat lower freeze-out temper- 
ature (AT ^ 20 MeV) at mid-rapidity as compared to protons. This difference in Tb is 
present throughout most of the phase space, only around target and projectile rapidities 
both measurements agree within errors. 

Similar results were reported by the EOS experiment [3] for Ni+Cu collisions at energies 
of about 2A GeV. Effective temperatures of 106 ± 5 MeV (A) and 142 ± 1 MeV (protons) 
have been derived from an exponential flt to the measured transverse mass spectra within 
the central rapidity range \y/ybeam\ ^ 0.25. 
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C. Production Yields 

The rapidity density distribution dN/dyc.m. is determined separately for each bin in ra- 
pidity by an analytic integration of Eg. [3] from = mo to ^ oc, using A and Tg 



parameters from the fit. One obtains |33l |: 



where mo is the particle's rest mass. Using this method, contributions to the integral from 
outside the acceptance are included. 

The resulting dN/dyc.m. distributions are presented in Fig. [71 The meaning of the different 
symbols is the same as in Fig. [6l The distributions for and A particles are corrected for 
the correspondin g br anching ratios. The upper plot shows a comparison between and 



distributions 



24 



34l |. which are found to be very similar. The total production yield 
per central (350 mb) collision of 0.078±0.003 and the corresponding width of the dN/dyc.m. 
distribution of 0.493±0.016 are obtained from a Gaussian fit to the data. The combined 
distribution has a width of 0.433±0.011 and a total production yield of 0.075±0.002. 

The dN/dyc.m. distribution of A particles is displayed in the lower panel of Fig. Hand is 
compared to the one of protons (down-scaled in order to have the same integral as the A 
data). One observes a pronounced longitudinal spread of the proton distribution which could 



be an indication for a large de, 
been found for other systems 



gree of transparency in the Ni + Ni system, which has already 



35I . 3^ . The A rapidity distribution is much more compact as 



expected for particles produced in the fireball. Such trend is also observed at AGS energies 
(2l| . The total A production yield per central (350 mb) collision is 0.137±0.005 and the 
width of the distribution is 0.386±0.009. 

The dashed curves in both plots denote a Siemens- Rasmussen distribution jsol (normal- 
ized to the respective total yield) for a temperature of T =90 MeV and radial expansion 
velocity (3 = 0.3 as found in [10]. While the kaons fully agree with the Siemens- Rasmussen 
distribution, a slight difference is visible for the A in the mid-rapidity region. This may 
indicate a non-isotropic emission of the A in the center-of-mass system. 
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D. Evaluation of the Systematic Errors 

In order to estimate the systematic error on the width of the rapidity density distributions, 
on the effective temperatures and on total production yields we apply different sets of cuts 
by varying the selection criteria of the relevant quantities discussed in section [IVl These 
selection criteria are determined so that it is possible to accept variations of counts in the 
invariant mass spectra of about ±50%. While the effective temperatures vary within 6-8% 
for both particles, the widths of the rapidity density distributions differ by about 6% (K^) 
and 12% (A). The total production yields change by about 10% (K^) and 18% (A). The 
systematic errors are summarized in Tab. [TTl 



systematic 
error of 




A 




+4 


+9 


slopes [MeV] 


-7 


-7 




+0.024 


+0.047 


dN/dy width 


-0.031 


-0.031 




+0.007 


+0.009 


total yield 


-0.008 


-0.025 



TABLE II: Systematic errors on the and A results. 

rn 

Furthermore, an independent analysis based on a multi-layered neural network [37], also 
has been used to identify A {32!, 38]. Within errors the results of yields and inverse slope 



parameters show a good agreement with the ones obtained from the standard procedure. 



E. Strangeness Balance 

In heavy-ion collisions, strange particles are produced via the strong interaction which 
conserves the total amount of strangeness in the system. At beam energies of about 2A GeV 
the bulk part of the produced strange particles comprises the A and S baryons and the kaons 
(X^, K^^ K^). Heavier strange resonances or multi-strange particles have [NN) production 
thresholds which lie significantly above this energy and thus are strongly suppressed. There- 
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fore, their contribution to the overaU strangeness production is neghgible and wiU not be 
considered in the foUowing. 

Due to the conservation of strangeness, the balanced sum of the production probabihties 
of the above mentioned particles then has to cancel: 

= S+ + S- + S V A V A^" + ^ - + . (6) 

Here the symbols of the particles denote the respective strangeness production probabilities. 
As mentioned earlier, the S° can not be disentangled from the A° due to its decay properties 
and to the detector itself. Thus one always measures the (A° + S^) yield. In the case of 
neutral kaons, one measures the decay of into two charged pions. After acceptance and 
efficiency corrections one obtains the yield of {K^ + K^) which is equal to {K^ + K^). 







stat. 


syst. 


particle 


yield 


error 


error 


K+ [24,^ 


0.075 


± 0.002 




K- [34] 


0.003 


± 0.001 






0.078 


± 0.003 


+0.007 
-0.008 


A + EO 


0.137 


± 0.005 


+0.009 
-0.025 



TABLE III: Measured production yields of strange particles in central Ni+Ni collisions (350 mb). 



Using the known yields of X+ [24, [s^ , y , X° + and A^ + in Ni+Ni collisions 



at 1.93^4 GeV beam energy given in Tab. Hill one can now determine the missing part of the 
bulk strangeness, i.e. the charged strange baryons S+ + S~. These baryons always decay to 
final states including neutral particles and thus can not be detected with the FOPI setup. 
In order to estimate the production yield of charged S baryons we use the following relation: 

S+ + S- = X+ + {K^ + K^)- + A^)-3'K- (7) 

under the assumption that the number of produced is approximately equal to the number 
of produced K~ . The result is a yield of 0.007 ± 0.008 (stat.) Iq^oit (syst.) particles per 
event (cf. section [V Dl for the determination of the systematic errors). Thus there is very 
little room for the production of charged (and due to isospin symmetry also neutral) S 
baryons. 



14 



This result can be understood partially from the fact that due to the mass difference 
between and S^'^ of about 80 MeV, the thresholds for the direct production in A^A^ 
collisions are different. While about 1.6A GeV beam energy is needed to produce a (A°, 
pair, the creation of a (S^'^, K^^^) pair requires energies higher than 1.8 A GeV. At a beam 
energy of 1.93yl GeV the A° therefore is produced at an excess energy (in the center-of-mass 
system) of about 125 MeV whereas only 50 MeV are available for the S baryons. Thus, one 
already expects a certain suppression of the with respect to the A°. 

Experimental studies of p-p collisions by the COSY- 11 collaboration [S] have shown 
that at (center-of-mass) excess energies of a few MeV, the direct production of the is 
suppressed by about a factor 30 relative to the A°. In another experiment by COSY-11, 



the energy dependence of the S+ cross section was investigated |40|] and was found to be 
enhanced by about a factor 10 relative to the A° results. Both, the suppression and the 
S+ enhancement seem to be a result of strong final state interactions between the outgoing 
hyperon and the remaining nucleon. 

Evaluating the results of COSY-11 for our situation of different excess energies for A° 
and S production, we find a suppression factor of about 30 for the S° with respect to the A° 
and an enhancement factor of three for the While the suppression is fully consistent 
with our results from the strangeness balance, the S+ enhancement is clearly incompatible. 
Transport and statistical model approaches [411] calculate a (S+ + S~)/ (A° + S°) ratio of the 
order of 0.35, while a value of 0.05 is found from our data. A possible explanation could be 
strangeness exchange occurring in collisions between S+ and neutrons with a in the final 
state. This channel would then increase the A° yield and produce the observed behavior of 
the strangeness balance. 

F. Centrality Dependence 

The yield per event of (full circles) and A (full squares) as a function of the centrality 
(i.e. impact parameter 6, cf. Tab. [T]) of the reaction is presented in Fig. [3 As expected, the 
yield per event increases strongly with increasing centrality, for both particles (top panel). 
One observes an almost linear dependence of the and A production yields on the impact 
)arameter. This has been already observed at BEVALAC in Ni + Cu reactions at 2 A GeV 



18[ |. Despite the qualitative agreement of both measurements for the A, our data show an 
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increase of the production yield from the most peripheral to the most central bin which is 
about 20% higher than that from the EOS experiment. The observed discrepancy lies within 
the systematic errors of the yields of both experiments. Due to missing errors for the mean 
impact parameters of the EOS centrality bins, the effects of multiplicity fluctuations and 
of the poor resolution of low impact parameters on the data are unclear. A more speciflc 
discussion of the observed discrepancy therefore is not possible. 

The total yield per event and per participating nucleon [P/ Apart) is plotted as a function 
of Apart (cf. Tab. [T]) in the bottom panel of Fig. El For both, and A, it exhibits a similar 
rise with Apart) which is of the order of 50% from the most peripheral (MUL4) to the most 
central (MULl) bin. The dependence of P/ Apart oii the number of participants can be 
parameterized by 

o,{A,,rtt-'^ (8) 



^part 

with = 1.20 ± 0.25 and = 1.34 ± 0.16 for and A, respectively. 

Within errors the data agree with results on the centrality dependence of the pro- 
duction yields measured by the KaoS collaboration for the same system [42] (star symbols). 
Note that for comparison, the KaoS data measured at 1.5 A GeV have been scaled by a fac- 
tor of 3.80 since it was shown that the increase of the production yield with increasing 
beam energy {E^^eam) is proportional to E^fa^''^ 13] • Similar trends were also evidenced for 
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and K~ measured at lower beam energies in Au + Au and Ni + Ni collisions 
The variation of the production yield per Apart as a function of Apart can be interpreted 
as a signature of multiple collisions which are dominantly contributing to the production 
of strange particles near threshold Q, [45I. Moreover, the similar trends evidenced for 
(K^) and A allow to conflrm that these particles are produced according to the same creation 
mechanisms: B B B Y K and n B Y where B stands for a nucleon or a A 
(A^*) resonance and Y for a A or a S hyperon, in the 1-2A GeV energy range jl]. 

Systematically higher values of a are expected when excluding the most central data point 
from the flt using Eq. El For the one flnds a^o = 1.38±0.54 and < = 1.84±0.35 for the 
A, respectively. The diflFerence between these results and the ones from the full flt mentioned 
above could imply a saturation of the strangeness production in very central collisions. The 
effects of the multi-step processes could be partially canceled by the reshuflling of strange 
particles by absorption in the large volume of the flreball. 
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VI. CONCLUSIONS 

We have reported on new results from the FOPI experiment concerning the production of 
neutral strange particles, and A, in Ni + Ni collisions at 1.93A GeV. A large statistics of 
high quality and A candidates, with a signal-to-background ratio of the order of one to 
two, have been reconstructed. Phase-space distributions have been measured down to low 
transverse momenta and in the full rapidity range. Differential yield distributions and inverse 
slope parameter distributions have been found similar for and K^^ confirming that their 
production mechanisms are identical and that there is no measurable effect of the isospin 
difference between the two particles. The rapidity distributions of protons and A exhibit 
different trends due to an incomplete stopping in the studied system. Regarding the overall 
strangeness production, the results request a strong suppression of the charged S hyperon 
production which is not seen in transport or thermal model calculations. The mechanism for 
this has yet to be found. For the first time at SIS energies, the total production yield of 
and A per participant is studied as a function of the centrality of the reaction. The increase 
of the yield with rising number of participants, similar for both and A, indicates that 
multiple collisions contribute to the production of neutral strange particles near threshold. 

All observed trends confirm that K^{K^) and A production mechanisms are strongly 
correlated. Detailed comparisons to transport models have to be done in order to check 
the consistency of the data with the hypothesis of modified in-medium properties of strange 
hadrons. 
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FIG. 1: (Color online) Schematic view of the and A reconstruction in the x — y plane (top) 
and in the r — z plane (bottom). The cut quantities used for the selection of the candidates are 
illustrated in the plots. 
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FIG. 2: (Color online) Invariant mass spectra of tt tt^ (top) and tt (bottom) pairs. The solid 
lines denote the combinatorics, the shaded areas represent the scaled mixed-event background. 




FIG. 3: (Color online) Phase space distributions of the reconstructed (top) and A (bottom), 
in the plane pt "^s. ycm.- The solid lines indicate the polar angular coverage of the CDC. 
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FIG. 4: (Color online) Reconstruction efficiencies for the (open circles) and the A (full squares) 
as a function of the transverse mass rut — tuq given for different bins of rapidity yc.m.- The dashed 
and solid lines denote fits to the efficiencies of and A, respectively. 
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FIG. 5: (Color online) Transverse mass spectra of reconstructed (left) and A (right). The 
spectra are plotted for ten rapidity bins ranging from yc.m. — -1-0 (lowermost spectrum, scaled by 
10^) to Hem. — 0.0 (uppermost spectrum, scaled by 10^). The data are corrected for the branching 
ratios of the respective decay channels and the unmeasured contribution. 
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FIG. 6: (Color online) Rapidity dependency of the inverse slope parameter for (full circles, 
upper plot) and A (full squares, lower plot). The open symbols are data points reflected with 
respect to mid-rapidity. The star symbols in the upper plot denote previously measured [2^ 
and the star symbols in the lower plot denote protons from our experiment. The lines represent a 
fit assuming an isotropic thermalized source (see text for details). 
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FIG. 7: (Color online) Rapidity density distributions for (full circles, upper plot) and A (full 
squares, lower plot). The open symbols are data points reflected with respect to mid-rapidity. The 
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34( 1 (upper plot) and 



(black and green) star symbols correspond to previously measured 
to protons from this experiment (lower plot). The lines represent a Gaussian fit from which the 
total production yield is extracted (see text for details). 
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FIG. 8: (Color online) Upper panel: Yield per event of A (full squares) and (full circles) as a 
"unction of the impact parameter h. The open squares indicate A data from the EOS experiment 
18|. Lower panel: Yield per event of A and per participating nucleon as function of Apart- The 

open stars denote experimental results for from the KaoS experiment 



